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Traffic behavior in the Pierce loop for data transmission is studied 
under assumptions of heavy loading. A deterministic mathematical model 
for describing traffic flows is developed and analyzed. The mathematical 
problem is of a linear complementarity type which has not been dealt 
with in the literature of mathematical programming. An effective procedure, 
the load-and-shift algorithm, for determining traffic flows is proposed. 
The procedure yields all feasible solutions for traffic flows and reveals 
the possibility of stations grouping into dominating classes and preventing 
other stations from using the system. This property, which can be eliminated 
by exercising appropriate control, also may affect the stochastic behavior 
of the system when heavy traffic conditions do not prevail and therefore 
deserves careful investigation. The paper includes two numerical examples 
illustrating use of the load-and-shift algorithm and numerical results 
from a simulation showing some of the effects of dominating classes when 
heavy traffic conditions do not prevail. 

I. INTRODUCTION 

The concept of a loop network for data transmission has been proposed 
recently by J. R. Pierce. 1 In such a network the stations are connected 
to a closed loop main line on which one-way traffic is allowed. A message 
to be delivered from one station to another is arranged, at the sending 
station, into standard packets each carrying the address of the receiving 
station. These packets are then delivered onto the main line, one at a 
time, where they flow around in the allowed traffic direction. The address 
of each packet is checked at each station on the way until it reaches 
the receiving station where it is removed from the main line. Traffic 
on the main line cannot be delayed; therefore, a station can deliver 
a packet onto the main circular line only when permitted by the existence 
of a gap in traffic or when receiving a packet from the main line. Principal 
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features of the system may be explained with the aid of the four- 
station network shown schematically in Fig. 1. 

The four-armed structure revolves around the central axis and stops 
briefly every time the four packet-carrying compartments at the ends 
of the arms are aligned with the four stations. During such a stop 
each station is able to check the content of the aligned compartment. 
If the compartment is empty, the station can load it with a packet. 
If there is a packet in the compartment, it will not be removed unless 
it is addressed to the said station, in which case the station is permitted 
to load the compartment again after unloading it. A Pierce loop can 
be represented by the mechanical analog structure shown in Fig. 1; 
however, the number of revolving arms in the structure is not neces- 
sarily equal to the number of stations in the loop. Rather, the number 
of arms is determined by the "loop time" of the system (the time needed 
for a bit to complete one round on the loop). Significance of the loop 
time is discussed in more detail in Section V. 

In the Pierce loop it is currently assumed that no outside control 
is applied and each station strives to send its messages at the earliest. 
Therefore, a station will never miss the opportunity to load a compart- 
ment unless there are no messages waiting for delivery at the station. 

This paper presents a study of the flow characteristics of such a 
system in heavy traffic, i.e., when the system is not able to deliver all 
messages, and infinite queues build up at some stations. In Section II 
the notion of "stable solutions" for the traffic flows is introduced and 




Fig. 1 — A schematic description of a four-station Pierce loop. 
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formulated mathematically. The case of a totally saturated system is 
analyzed in Section III where the basis is laid for the mathematical 
development of the load-and-shift algorithm given in Section IV. A 
detailed description of a Pierce loop for data transmission is given in 
Section V where the alternating priorities effect due to dominating 
classes of stations is studied by simulation. This section also discusses 
the important aspect of the order of stations in the loop. Readers not 
interested in the mathematical elaborations may skip Section IV. 

II. MATHEMATICAL FORMULATION 

We assume that the flow direction on the main circular line is counter- 
clockwise. There are n stations, connected to the main line, numbered 
from 1 to n in counterclockwise increasing order. The segment of main 
line between the ith and the (i + l)th station, i = 1, 2, • ■ ■ , n — 1, 
is called the ith branch. Similarly, the nth branch is the segment be- 
tween station n and station 1. Let p i{ be the proportion of flow (packets) 
emerging from station i and destined for station j. 

Ep«-i. (i) 

1 = 1 

The n X n square matrix P = [pa] possesses all the properties of a 
stochastic matrix. Note, however, that the elements of P are not 
necessarily probabilities. 

The demand at station i is given by X, which is the average amount 
of flow (packets per time unit) generated at the station. The capacity 
of each branch equals 1, that is, each branch is capable of carrying 
a maximum flow of one unit. In the schematic description given in 
Fig. 1 assume that one full revolution takes four time units (in general 
it will take the number of time units equal to the number of arms) 
The capacity of each branch will then be one packet per time unit. 

The average flow emerging from station i will be denoted by x { 
and the average flow in branch i by p, . Clearly 

S Xi ^ X, : , i = 1,2, ■•• , n (2) 

and 

S Pi £ 1, i = 1, 2, • • • , n. (3) 

Let a,, be the proportion of flow emerging from station j and flowing 
through branch i, then 

a,, ■ = Z Pi* i ( 4 ) 

ktSij 
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where 



(i+ 1, » + 2, ■■■ , j) 



if j > i, 



8 ti = -T ' *" ' "' "' "" ' ' "' (5) 

[(i + 1,»+ 2, ••• ,n, 1,2, ••• ,i) if j ^ i. 

Note that a t _i,— a,, = p,-,- ^ for all j 5^ i where i— 1 is defined to equal 
n when i = 1. For j = z we have a,, = 1 and a.-n = p,, ^ 1. In most 
reasonable applications p,< = 0. 

The average flow in branch i may now be expressed as a linear func- 
tion of X — (xi , x 2 , ■ • • , x„). 



>■• = Zl a u x i = 1 • 



(6) 



Every X which is a feasible solution for the average flows must satisfy 
relations (2) and (6). The set of all feasible solutions is therefore con- 
tained in a convex polyhedral set. A central control could select a 
particular solution from this set to suit a given objective. 

In the circular network suggested by Pierce, however, there is no 
central control. Rather each station is striving to maximize its own 
flow onto the main line. For this case we define a stable solution as 
a solution from which the system will not depart without outside 
intervention. Suppose then that X* is a stable solution. Clearly X* 
must satisfy relations (2) and (6) and the additional condition that if 
p, = 2^"_i a ti x*- < 1 then x* t = A, . To show that this is a necessary 
condition assume that p, < 1 and x* t < X, . However, station i strives 
to maximize its flow and can increase it as long as x* { < X, and p, < 1. 
Therefore, it is not possible that p, < 1 and x* { < X, . This additional 
condition is not generally sufficient for assuring that x* is a stable 
solution. 

Insufficiency is best demonstrated by a simple numerical example. 
Suppose n = 4, Xi = X 2 = X 3 = X 4 = 2 and 



P = 



10 
10 
1 











It is easy to verify that X° = (0, 0, 1, 1) and X + = (1, 1, 0, 0) both 
satisfy relations (2), (6), and the additional necessary condition as 
does X = aX° + (1 — a)X + for0 ^ a ^ 1. For simplicity we assume 
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that the flows generated at the stations are deterministic in nature 
and each station generates exactly two packets each time unit. 

Suppose now that the system (see Fig. 1) has one rotating arm only. 
In such a case there are only two stable solutions, namely X°, (a = 1) 
and X* , (a = 0). If the system has two rotating arms, an additional 
stable solution, a = \, is added. If the system has k arms, there exist 
k + 1 stable solutions, a = m/k, m = 0, 1, • ■ • , k. The system shall 
settle for the stable solution a = m/k if at time zero m compartments 
contain packets from stations 3 and 4 and k — m compartments contain 
packets from stations 1 and 2. The necessary condition for stability 
is also a sufficient condition if the flows generated by the stations 
are continuous (as in the case of nonmixable fluids or small particles 
such as vehicles and a loop consisting of a pipe or a road). In such a 
case the set of feasible solutions is identical to the convex polyhedral 
set given by relations (2) and (6). 

We wish to find all solutions which satisfy relations (2), (6), and 
the necessary stability condition (these will be all the stable solutions 
in the case of a system with continuous flows). Our problem can be 
redefined as one of finding all feasible solutions to the following set of 
equations: 

x t + Ui = X, , 
I>„s, +Zi = 1, i = 1,2, ••■ ,n. (7) 

m,z, = 0, .r, ^ 0, ^ ^ 0, z, ^ 0. 

This form resembles a linear complementarity problem, 2 where w, and 
2, are slack variables. A feasible solution to the set of equations (7) 
is a basic feasible solution to the set of equations (2) and (G) since at 
least n of the variables must equal zero. 

III. COMPLETE SATURATION 

For sufficiently large values of X, (for example X, = 1, i = 1,2, • • • , n), 
all the branches are saturated and p, = \, i = 1, 2, • • • , n. The set 
of equations (7) takes the form 

i>„,r, = 1, i= 1,2, ••• ,n. (S) 

/ = i 

This can be shown to be equivalent to the set of (n + 1) linear equations. 
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n 

Xi — 23 Pa x > = 0, * = 1, 2, • • • , n, (9a) 

E a;, d, = 6. (9b) 

j=i 

The ith equation in (9a) is obtainable by subtracting the ith. equation 
of set (8) from the (i + l)th equation of set (8). Equation (9b) may 
be selected as any linear combination of equations (8). Note that one 
equation in (9a) is redundent since 23? -i Va = 1 f° r all *'• 

Remark: The physical interpretation of equations (9) is that when 
in complete saturation a station is able to deliver only when receiving. 
The flow emerging from a given station must equal the flow entering 
the station from the main line. This equilibrium relation is expressed 
by the n equations of set (9a). The (n + l)th equation, (9b), expresses 
the capacity limitation of the branches. If we select d, = 2^v=i a a 
we have 

23 Xi dj = n, (10) 

»=i 

where dj is the average number of branches (distance) traveled by a 
packet emerging from station j. Equation (10) states that the average 
work (in terms of packets times distance) demanded from the system 
per time unit must equal n, since all branches are saturated and each 
traverses one packet per time unit. In matrix form we have 

X T P = X T , 
Z «»• *i - b> (11) 

1 = 1 

where X T is the transpose of X, (note that all vectors are defined to 
be column vectors). 

P is a stochastic matrix and, therefore, the problem represented 
by equations (11) strongly resembles one of determining the steady 
state probabilities of a finite state space Markov chain. The difference 
is that in the Markov chain problem dj = b for all ; while in our problem 
this is not necessarily so. In the following we shall make use of this 
resemblance. 

Definitions: 

(i) Sj shall be used to abbreviate "station i." 
(ii) Sj is said to be accessible from S { , S { — > Sj , if there exists a 
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sequence of elements of P such that P ik ,P k ,k, • • • Pkr-,k r Pk r j > 0. 

In such a case it is also customary to say that *S, leads to Sj . 

(Hi) If Si —* Sj and S, — > <S, both stations are said to communicate 

(Si «-> Sj). Clearly if S t <-> S, and S k <-» S, then & <-» S* . 

(iv) Let C(i) = («S, : S f <-> S t }. Clearly if S t t C(i) then C(j) = C(i). 

(v) A nonempty class of stations, C, is called a communicating 

class if for some station S, t C, C = C(i). It follows that two 

communicating classes are either identical or disjoint. 

(vi) A communicating class is closed if no station outside the class 

is accessible from a station in the class. 
(vii) A closed communicating class shall be called a dominating class 
or a class of dominating stations. A station not belonging to any 
dominating class is a dominated station. 

From the theory of Markov chains we know that there exists at 
least one class of dominating stations in a given Pierce loop. However, 
it is possible that there will be no dominated stations. In such a case 
all the stations are dominating and may form into one or more dominat- 
ing classes. It may be shown that the number of dominating classes 
is one if and only if there exists a station accessible from all other 
stations in the loop. In the case that dominated stations exist each 
must lead to at least one dominating station. 

Each dominating class of stations is represented by a principal sub- 
matrix of P (a dominating submatrix). This submatrix is obtainable 
by deleting all rows and columns of P corresponding to stations not 
belonging to the particular dominating class. Similarly, all dominated 
stations may be represented by one submatrix of P. 

Theorem 1. Let B be a k X k submatrix of P representing a dominating 
class of stations, then all vectors Y satisfying the equation, 

Y T B = Y T (11a) 

form a linear space of one dimension. Furthermore, all the elements of Y 
must have the same sign, i.e., either all positive, or all negative, or all zero. 

Proof. Theory of finite Markov chains (e.g., Kemeny and Snell 3 ). 

Corollary. Assume that the dominating class represented by B is C = 
{ Si, , St, , • • • , Si b ) then there exists a unique vector Y* = (?/ f * , y* , • ■ ■ , 
y* ik ) satisfying equation (11a) and the scaling equation 

£ Vt di, = b. (12) 



2528 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1971 

Clearly Y* > and any solution to equation (11a) may be obtained 
by multiplying Y* by some real number. We define a vector X = 
(xi , x a , • • • , x n ) where x { = if 8 t 4 C, and z, = y* if i = i t . The 
vector X is called a dominating solution to equations (11) corresponding 
to the matrix B and the class C. 

Theorem 2. Let Q be a k X k submatrix of P representing dominated 
stations, then the only solution to Y T Q = Y T is Y = 0. 

Proof. Theory of finite Markov chains. 

Theorem 3. Suppose that the matrix P contains exactly m dominating 
submatrices B x , B 2 , ■ • • , B m representing dominating classes d , C 2 , 
• • • , C m and suppose that X, is the dominating solution corresponding 
to B t , i = 1, 2, ■ • • , m. X is a solution to equations (8) if and only if 



and 



X = 22 &\Xi , 



2>. = 1. (13) 



ivhere a, , i = 1,2, • • • , m are real numbers. 

Proof. This is an immediate result of Theorems 1 and 2. 

We conclude that in the case of complete saturation, (p, = 1, i = 
1,2, • • • , n), there always exists a nonnegative solution to our problem. 
If there exists only one dominating class of stations in the loop, there 
exists a unique solution to equations (8). If there is more than one 
dominating class in the loop, then there exist infinitely many non- 
negative solutions obtainable as convex combinations of the dominating 
solutions. 

IV. THE LOAD-AND-SHIFT ALGORITHM 

Returning to the more general case, we describe in this section the 
load-and-shift procedure for solving the set of equations (7) when 
Pi S 1, i = 1,2, • • ■ ,n. The algorithm is based on the results obtained 
in the foregoing analysis. Assuming that the capacity of the branches, 
e, may be varied between and 1, we start with e = and increase 
it until e=l. While so doing we simultaneously load the system and 
obtain the feasible solutions for any given value of ^ e ^ 1. 

We shall start by outlining the procedure for finding just one solution. 
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Start: r: = 0, e: = 0, P (1) : = P, X'": = X, . 

Step 1: Increase r by 1. In P (r) find a square submatrix B <r) repre- 
senting a dominating class of stations C r , C r = {S il , 
£,,,••• ,S it }S 

Step 2: Find the unique positive vector F (r> = [y\? , y£ } , • • • , y^] 
satisfying the set of equations 

(Y (r) ) T B (r) = (Y {r Y, 

To determine the value of d im we select an i such that S t t C (r) 
and let d im = a,, lm . 

Enlarge F (r) to the form X (r) = (x{ r) , x { 2 r) , ■•■ , x l n r) ), 
where x< r> = if S,- * C (r> and x) T) = yV* if j = i m . 
Step 3 (Load): Find a number A <r) such that 

where the minimization is over all i such that S { t C T) . 

If A <r> ^ 1 - e, set A (r) = 1 - e, set iV = r, and go to "Last 
Step." Otherwise increase e by A (r) and continue. 
Step 4 (Shift): If r = w, set 2V = n and go to "Last Step." Otherwise 
construct the (n — r + 1) X (n — r + 1) square matrix P (r+1) 
by adding the ;'th column of P (r) to its (j + l)th column (if j 
is the last column, it is added to the first one) and then deleting 
the jth column and the jth row. Xj ,+1) : = X, <r> - A (r, x, (r) , 
i = 1, 2, • • ■ , n. End of rth iteration. Go back to "Step 1." 

Last Step: The solution is 

X = £ A (r) X (r) . (14) 

r = l 

STOP. 

Theorem 4: The procedure described above will always yield a vector 
^ X ^ X m of ??*os£ ?i iterations. 

Proof: All possible matrices P (r) are stochastic. Therefore there always 
exists at least one dominating submatrix of P <r> denoted by B ir) . 
Since X (r) ^ and A (r) ^ then X = ^ ml A (r) X (r) ^ 0. From the 



A procedure for determining B lr) is described in the Appendix. 
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algorithm we have X = £*i A (r, Z (r) + \ lm = X + X (Ar> . The algo- 
rithm also ensures that X (r) ^ 0, r = 1,2, • • • , N, and therefore X ^ X. 
Since P is n X « the number of iterations cannot exceed n. 

Theorem 5 {Existence): X obtained by the Load-and-Shift procedure 
is a feasible solution to equations (7). 

Proof: We enlarge F (r) by adding zero elements corresponding to 
columns of P (r) not included in B iT) . The enlarged vector, denoted 
by Y olr) is a dominating solution to 

Y T P (r) = Y T 
Z to * - 1, (15) 

jtBr 

where J? P is the set of indices of columns included in P (r) . These equa- 
tions are equivalent [see equations (8) and (9)] to 

E aifvt : = 1, iefl r . (16) 

Jtflr 

It is easily verified, by the use of equation (4), that 

<*? = a,,- . (17) 

The vector X lT) is obtained by adding to F 0(r) zero elements correspond- 
ing to columns of P not included in P (r) . Therefore 

A ,Tl Z^! rl =A (r ', for all i t R r . (18) 

From the definition of a,, , equations (4) and (5), we know that 
Q>i-u ^ 0-a for all j except j = i, (where i — 1 = n if i = 1). If i $ R r 
then x- r) = 0. It follows that 

A (r) Eo iA w ^ A (r) , for all »>B r . (19) 

Summing equations (18) and (19) with respect to r and then substitut- 
ing equation (14) we obtain (for N < n) 

n N 

Z a«xt = Z A (r) = 1, ^ * * Rr for all r = 1, 2, ■ ■ • , N, 

1=1 r=l 

n AT 

Z a>u%i ^ Za (,, = 1, otherwise. (20) 

j=l r=l 

From the algorithm we know that if i t R r for all r = 1,2, • • • , N then 
Xi ^ X, . Otherwise x { = X, . This completes the proof for N < n. 
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In the case N = n it is possible that £jl, A <r) < 1 (nonheavy traffic). 
The theorem still holds since x { = X, , * = 1,2, • • • , n. In the following 
we illustrate the use of the algorithm in solving for the flows in a 5- 
station loop: 



P = 
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The procedure terminated in three iterations yielding a solution: 

Xi = 0.4646 

x 2 = 0.4000 = X 2 

x 2 = 0.5000 = X 3 

x 4 = 0.2333 

x s = 0.0508 

Actually the algorithm finds the value of one variable in each iteration. 
When all variables equal to their respective X,'s have been determined, 
the algorithm finds in one iteration the values of all the remaining 
variables. It is therefore advisable to test the possible solution X = X 
beforehand. 

We note that in the example there is a single dominating class of 
stations in each iteration (i.e., each matrix P (r> has only one dominat- 
ing submatrix B (r) ). It will be shown later that if P has only one dominat- 
ing submatrix then there exists a unique solution to equations (7). 
In the general case, however, P (r) may have more than one dominating 
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submatrix. This may give rise to the existence of infinitely many solu- 
tions. We would now like to improve the algorithm in order to be 
able to determine all feasible solutions. 

Theorem 6: If P (r) has exactly one dominating submatrix, say B lr) , 
then P <r+I) will also have exactly one dominating submatrix. 

Proof: Suppose i is such that S t c C <r) , then S t — » S t for all j. We 
apply the shift operation from i to k, (adding column i to column k 
and then deleting row and column i), thus creating P (r+1) . Clearly 
now Si — * S k for all j. However, if there exists a station accessible 
from all other stations there exists exactly one dominating class in 
the loop. Furthermore C (r+1) = C(k). 

Corollary: If P has one dominating submatrix so will P (r) , r — 1, 2, 
■•■ ,N. 

Theorem 7: If P (r) has m > 1 dominating submatrices, P <r+1) will have 
either (m — 1) or m dominating submatrices. 

Proof (outline): Let B[ T) , B 2 T) , ••• , Pl r) be the m dominating sub- 
matrices of P (r) and let C[ T) , C 2 \ ■ ■ ■ , Cl r) be the corresponding 
classes of dominating stations. Q (r) is the submatrix representing 
dominated stations. We apply the shift operation from i to k to obtain 
P (r+1> . Without loss of generality we assume that £, e C\ T) . For S k 
one of three alternatives must be true: 

(i) S k e C[ r) . In this case a dominating class will be formed, con- 
taining some or all the remaining stations of C[ T) . Note that S k must 
be in the newly formed dominating class since it is accessible from all 
remaining stations of C[ r) . Those stations of C[ T) which are not included 
in the newly formed dominating class turn into dominated stations. 
The matrix p (r+1) will then have m dominating submatrices, namely 
B l 2 r) , B ( 3 r) , - • • , B^ } and a newly formed one. 

(it) S k £ C\ r \ I 9* 1. In this case every station in C\ T) is accessible 
from any remaining station of C\ r) . Therefore all remaining stations 
of C\ T) become nondominating, and P (r+1> will have m — 1 dominating 
matrices, namely B 2 r) , B l 3 r) , • • • , Pl r) . 

(Hi) S k is nondominating. S h must lead to at least one dominating 
class. If it leads to any dominating class other than B[ T) then this 
case becomes the same as (ii). If S k leads to B\ r) only, then a new 
dominating class is formed. This class includes all stations accessible 
from S t . Remaining stations of B\ r) which are not accessible from S k 
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become nondominating. P (r+1) will have m dominating submatrices, 
namely B { 2 T) , B { 3 r> , • • • , B% ] ', and a newly formed one. 

This completes the outline of the proof. It is important to note that 
a dominating submatrix of P (r) , say B\ r) , will be a dominating sub- 
matrix of P (r+1> if S, ; 4 C\ r) . 

An Outline of the Complete Version oj the Load-and- Shift Algorithm 

In the procedure for finding just one solution, outlined in the pre- 
ceding, we apply one load and one shift operation in each iteration. 
In the complete version of the algorithm we need to apply several 
such operations in each iteration. To eliminate possible confusion the 
superscript denoting the iteration number will be placed at the upper 
left side. Thus, for example, <rl P is the stochastic matrix remaining 
by the beginning of the rth iteration. 

Definitions: Suppose (r) P has m r dominating submatrices, lr) Bi , 
lr) B 2 , ■■■ , {T) B mr , representing classes (r> d , <r) C 2 , • • • , (r) C mr . 
The matrix <r) P was obtained from U) P = P by executing some sequence 
of load and shift operations. In order to keep track of these operations 
we relate to tr) B { a set <r) E t containing labels. If the label w a t e ir) E { 
we know that <r) -B, cannot be obtained unless the appropriate shift 
operation is applied to a) Bj . Thus, when (r) /?, , i = 1, 2, ■■ ■ , m r , 
are given, the exact sequence of load and shift operations that led us 
to (r) P is known. 

The n dimensional vector (r) X denotes the remaining unutilized 
flows at the n stations of the loop. The m T dimensional vector <r) e = 
( (r) «i , <r> «2 , ••• , <r) «m r ) describes the amount of branch capacity 
utilization. Thus <r) €, is the amount of branch capacity utilized by 
the sequence of load and shift operations resulting in <r) Z?, . 

For each (r) -B, it is possible to determine the dominating solution 
1t) X { and the quantity (r) A, in the manner described in steps 2 and 3 
of the procedure outlined previously. 

We say that ir) B { -* (r, fi y if when applying a shift operation on the 
appropriate column of (r) P, all remaining stations of <r) C, lead to 
stations of <r) C, . If (r) C, consists of a single station then (r) J5, — * <r) 5, 
if the shift is from this station to a station in (r) C,- or to a dominated 
station leading to <r) C, . The submatrices (r) -B, , i = 1, 2, • • • , m T , 
divide into three types: transient, terminal, and ring members. 

Ring: The set R = i <r, B„ , (r, -B„ , ••• , (r, B, t } is called a ring if 
it satisfies three conditions: 



2534 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1971 

(t) M B ix -> M B U -» > lr) B it -> M B U , 

(ii) Bi t , j = 1, 2, • • • , k, does not lead to any dominating 

submatrix outside of the same ring, 
(Hi) 1 - ZU ( Cr> «, - (r, A,,) > 0. 

Terminal: (i) If conditions i and ii above are satisfied while con- 
dition Hi is violated the matrices ir) B il , ir) B it , ■ • • , lT) B ik 
are called terminal matrices. 
(ii) If (r, A, ^ 1 - (r) e, then <r, 5, is terminal. 
(Hi) If lr) B { is terminal and (r) P,- -» <r) B, then <r) S, is 
terminal too. 

Transient: (r) S, is called transient if it does not belong to a ring 

and is not terminal. 

There is a strong similarity between the notion of a ring and the notion 
of a dominating class. To identify rings one may use essentially the 
same technique proposed in the Appendix for identification of domi- 
nating classes. 

It is important to observe that when applying the appropriate shift 
operations to all members of R, what remains of R will form one domi- 
nating submatrix and possibly some columns and rows corresponding 
to dominated stations will be left out. This property follows from 
Theorem 7. 

The following is an outline of the algorithm: 

Start: r: = 0, (1) P: - P, (1) X: = X, (1, e: = 0, (1) £, : = <f>. 

Step 1 : r: = r + 1. 

Find all the dominating submatrices of <r> P denoted by (r) Si , 
ir) B 2 , • • ■ , lr) B mr . Calculate the numerical values of (r) A, 
and (r> X, , corresponding to (r) B, , <r> X, and (r, e, ,i = 1,2, • • • , 
m r . Find all the rings of <r, P denoted by <r, #, , ir) R 2 , ■•• , 
{r) R kr . Identify the terminal and transient dominating sub- 
matrices of <r) P. If all dominating submatrices of <r) P are 
terminal go to LAST STEP. If the number of columns of 
(r) P equals the number of its dominating submatrices who 
form into one ring, go to LAST STEP 1. Otherwise <r) A, : = 
for all (T) B t not belonging to a ring. 

Step 2: Execute the appropriate shift operations to all (r) B { belonging 
to rings and obtain (r+1 P. Note that in (r+1) p there will be 
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one dominating submatrix corresponding to each ring of (r) P 
and, in addition, all terminal and transient dominating sub- 
matrices of <r) P will be dominating in <r+1) P. 

Step 3: 

" + » X : = ( "x - £ ">A, (r) X< . 



<,:= E ( <r, 6,- + (r, A ; ), <- 1,2, ••• , k r . 



(r + l) 

e, . = 

"+» € . : = <" 6j . where (r + ,) 5, = lr) B t , i = fc r + 1, k r + 2, • • • , 

m r+1 . 

(r+,) £, : = U C% U (r, a,), t = 1, 2, • • • , fc r . 

I; :<'>B,t<"R,l 

i**»B { : = ">#. where (r+,) S, = (r, £ y , i - k r + 1, fc r + 2, • • • , 
m r+1 . GO TO STEP 1. 

Last Step 1 : The unique solution is X = X. STOP. 

Last Step: N: = r. 

All the solutions are given by 

X = E E <"«, M X t , (21) 

r=l i=l 

E E (r «. - 1, =g '"a, ^ (r, A. , <"«, = 0, r £ 2, 

r - 1 i : - 1 

unless a> a, = U, A, for all '*'«,• e (r> #,- . STOP. 

Note that it is advisable to test the possible solution X = X beforehand. 

Theorem 8: The load-and-shift procedure will yield vectors ^ X ^ X 
in af n?os£ ?! iterations. Each such vector, given by equation (21), is a 
feasible solution to set of equations (7). 

Proof: The proof is practically identical to the proofs of Theorems 
4 and 5. 

The use of the algorithm is illustrated by a 10-station numerical 
example. 
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It remains to show that the load-and-shift algorithm will yield all 
the feasible solutions to equations (7). 

Theorem 9: If X° = (x\ ,x\ , • • • , z°) is a feasible solution to equations (7), 
then 

X, ^ «?£ Min|Sp^/,X*}' (22) 

Proof: The left-side inequality of (1) is part of equations (7). For the 
right-side inequality, suppose first that x° < X, , then 

^2 a f jx1 = 1 and ^ a.-jyx" ^ 1. 

Taking the difference we obtain 

E (a,, - a^M ^ 0. (23) 

From the definition of a,, we have that 

a,, — a,_i,- = —pa , for all j ^ i, 
and 

a,, — tti-u = 1 — Pa ■ 
Therefore 

*5 fc t, P,.*? . (24) 

It follows then that x] = X, if ^" =J p,-,x° ^ X, and therefore 

X, ^ z° ^ Minji;?),,^ , X, 

Theorem 10: Let X° be a solution to equations (7), and assume that X > 0, 
then 

(i) If x° = 0, then z° = for all i such that S, -* £, . 
(tt) // x) > 0, then x° > for all i such that S,- -> 8 { . 

Proof: If &, — » *S, there exists a sequence p,-, fe , p* lltl , • • • , P* r -,.*r > 
p*,,, whose product is positive. Suppose x] > 0, then min (2^_, Pm*,^-! > 
X*,) ^ Vikxtf > 0- From relation (22) it follows that x n kl > and simi- 
larly x° k% > 0, • • • , x" > 0. This proves ii. The proof of i is immediate. 

Remark: The assumption X > is not restrictive. If X, = 0, we apply 
a shift operation to the ith column of P. 
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Theorem 11: Let B (l) be a dominating submatrix of P and assume, 
without loss of generality that C {i> = \S t , S a , • • • , jS*} is the class of 
dominating stations represented by 5 ( ". A ( " and X ll) are then uniquely 
defined. Suppose X° is a feasible solution to set of equations (7) then 

Either: (i) x° = ax\ l) , i - 1, 2, • ■ ■ , k, £ a £ A (,) . 

or: (ii) x" ^ A (1) x[ n , i = 1, 2, • • • , k, where the inequality 

is strict for at least one value of i. 
(Hi) If i is the case then x" = for all Sj i C (I) and leading 

to C (l) . // there exists such an S t and x] > 0, thenii is 

the case. 

Proof: Since B (1) is a dominating submatrix then ]£*_, p,-,- = 1 for 
i - 1, 2, • • • , k. 

(i) We assume that x° < \< for i = 1, 2, • ■ • , k and then sum the 
first k inequalities (24) to obtain 

i=l i = l i-fc+l i-l 

For this relation to hold it is necessary that x" ^2 k i=1 p,, = for all 
j ^ k + 1. If C (1> is not accessible from £, then _>_)*., p,-,- = 0, other- 
wise cither ^* =1 Pa > and therefore x° must equal zero, or S t — * S,„ , 
m ^ k + 1, and 5_X, p mi > and therefore x„ must equal zero. From 
Theorem 10 we know that in this event x] = 0. Concluding then that 
x" = for all (j: *S, i C ( " | and S t — > C (I) we obtain [in a manner similar 
to the one used in obtaining relation (24)] 

k 

x°i = _Cp,-,a£ , i = 1,2, ■ ■ ■ , k, 
i-i 



yielding 



i = 1, 2, • • • , k. 



Since ^ .r° < X, we have that ^ a < A (n . It also follows that if 
there exists £, i C (1) and leading to C ( " and x" > then there exists 
at least one value of i, i = 1,2, • • • , k, for which x° = X,. 

(ii) Suppose there exists a value of i, i = 1,2, • • • , k, such that 
xl ^ A ll) x\ 1} . From Theorem 10 it follows that x • > for all i = 
1, 2, • • • , k. We select a number a, A (1> > a > 0, such that 

Min {sS-ozH = x?_ -aa^ = 0. 

I ■= 1 , 2 , • • • , k 

Clearly if such an a exists then x° im < A U) .tV b and from relation (24) 
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we have 

x° im ^ £,p iim x1. (25) 

7=1 

On the other hand we have 

X° im ««£ = Ep,^' 1 . (26) 

Subtracting equation (26) from (25) yields 

o= tp^-^ru £ ?,<„*?. 
i-i ,- - k 1 1 

Since S in t C u> there exists a value of j, j 7^ i m and £,- t C (1) , such that 
p iim > 0. It follows then that x\ = txx^ for i - 1, 2, • • ■ , m. This 
is a contradiction, which completes our proof. 

Corollary: Let ll) Bi , ll) B 2 , • ■ ■ , a) B mi be all the dominating submatrices 
of P = (1) P, and let ll) C t , (1) C 2 , • • • , (1) C mi be the corresponding domi- 
nating classes. If X° is a solution to equations (7) then either 

X n = £ (,) a,. W X< , Z <% = 1, ^ (,) a,- ^ (I) A,- , (27) 

i=i i=i 

or there exists at least one class of dominating stations, say (l) Ci , such 
thatX ^ ft V% . 

Theorem 12: If X° is a feasible solution to equations (7) it is obtainable 
by the load-and-shift algorithm. 

Proof: If X° is given by equation (27), it is obviously obtainable after 
one iteration of the algorithm. Otherwise there exists a class of dominat- 
ing stations, say U) C, , such that X° ^ w ^ t t "X l . We execute a 
shift operation involving the submatrix (1I J3, and obtain a stochastic 
matrix P <2) (note that P (2) * (2) P). We let X 0l2) = X° - (1, A, (,, X,. 
and X <2) = X - '"A,"'!, . Clearly Z 0<2) is a feasible solution to equa- 
tions (7) when using as parameters P <2) and X <2) and replacing the 
right side by 1 — (1) A, . It is possible therefore to proceed with a 
sequence of load-and-shift operations until X° is obtained. Since the 
algorithm takes into account all possible sequences of load-and-shift 
operations, X° is contained in the set of solutions given by equation (21). 

Corollary {uniqueness): If P has exactly one dominating matrix then 
there exists a unique feasible solution to equations (7). Note, however, 
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that this is not a necessary condition for uniqueness. It is possible that P 
will have several dominating submatrices and the solution will be unique. 



v. DISCUSSION 



The heavy traffic assumption enables us to regard the system as a 
deterministic one. The analysis of the deterministic flows shows that 
stations tend to "band" into classes with the ability of dominating 
the system and preventing other stations from using it. 

This undesirable property, which can be eliminated by exercising 
appropriate control, also may affect the stochastic behavior of the 
system when heavy traffic conditions do not exist. One can imagine 
two classes of stations competing for domination of the system. Since 
traffic is not heavy, all the stations in the system are able to deliver 
their messages in finite time. Nevertheless, when one dominating class 
controls the system, it will prevent other stations from using the belt 
line until the queue at one of the stations belonging to this class becomes 
empty. At that moment the competing dominating class is able to 
take over and prevent other stations from using the belt line. This 
may result in a situation of alternating priorities (see Ref. 4) where, 
while one class is served, the queues at the competing class build up. 
While average queue sizes may not be strongly affected, the strong 
fluctuations in queue lengths may be undesirable. This possibility has 
been explored numerically by the use of a digital simulator. 5 The 
operating principles of the simulated system will be explained with 
the aid of Fig. 2. 

Each of the stations is represented by a B-box. A packet coming 
out of a station is first multiplexed on the line by the B-box, provided 
the line is free, and then is passed from B-box to B-box until its destina- 
tion is reached. At each B-box on the way the address of the packet 
is examined. At the particular B-box of destination the packet is taken 
off the line. The main function of the A-box, shown in Fig. 2, is syn- 
chronization of the loop. 

Assume that the packets are made of L bits each and the address 
is given in the first k bits of the packet. A time unit in this system is 
the time it takes to multiplex a bit on the main line by a B-box. Assume 
also that the traveling time from one B-box to the adjacent one is 
zero. Station i is allowed to start sending a packet at times mL + ik, 
in = 0, 1, 2, • • • , providing that the main line is free. The A-box is 
a buffer. Bits coming out of the B-box of station n accumulate in the 
A-box. Suppose at time k B-box 1 starts sending a packet. At time 



2544 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1971 



B' 



STATION 

4 




STATION 

2 



Fig. 2 — A schematic description of a loop with n stations. 

k + 1 the first bit reaches B-box 2 and is delayed there until time 2k 
when the whole address of the packet has been received. If the packet 
is addressed to station 2, it will be taken out and at the same time 
B-box 2 can start sending out its own packet. If the packet is not ad- 
dressed to station 2, it will be sent from B-box 2 to B-box 3 where the 
same process will take place starting at time 3k. Bits arriving at the 
A-box are buffered. At time L the A-box starts to send bits (at the 
same rate as a B-box) until a whole packet has seen sent. If the buffer 
is empty the A-box will wait another L time units and will start sending 
at time 2L. In general, the zth B-box checks its buffer at times mL + ik, 
m = 0, 1, 2, • • • , and if the buffer is empty it may start sending its 
own packet. If the buffer contains the address of station i, the B-box 
will remove the arriving packet and may, at the same time, send out 
its own packet. If the buffer contains an address different than station 
i, the B-box will pass on the arriving packet. 

In a similar manner the A-box checks its buffer at times mL, m = 
1, 2, ■ • • ,. If the buffer is not empty, it sends out L bits. If the buffer 
is empty, the A-box remains inoperative for the next L time units. 

The loop time is the time it takes a bit to complete one round of the 
loop, and is measured in multiples of L. In the single A-box loop de- 
scribed here the loop time is the smallest integer greater than or equal 
to nk/L where n is the number of B-boxes in the loop. Clearly there 
is a complete analogy between the loop described here and the one 
presented in Fig. 1 if the number of revolving arms is taken as equal 
to the loop time. 
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It is important to note that the simulated system described here 
represents only one conceptual way in which the loop may be operated. 
It is possible, for example, that the packets will move from one B-box 
to another in one block rather than bit by bit. This is equivalent to 
placing an A-box between any two B-boxes (n rotating arms). 

We have used the digital simulator to examine the queuing char- 
acteristics of several small systems. The main purpose of the simula- 
tion was to study the effects of dominating classes in nonheavy (non- 
saturated) traffic situations. Numerical results are presented for an 
8-station loop with two dominating classes d = (aS x , S 2 , S 3 , S t ) 
and C 2 = (S 5 , S 6 , S 7 , S s ). The loop time for this system was selected 
to equal 1 (one rotating arm) and the system was simulated for three 
different expected main line loads (utilization). The P matrix and 
average queue sizes (in packets) at the stations are shown below. 

Simulated Example: 8- Station Loop (Two Dominating Classes) 



P = 



S, 

S 2 
S* 
8 4 

s 5 
s- 

s 7 
s 



Si o 2 S 3 Si 

1/3 1/3 1/3 

1/3 1/3 1/3 

1/3 1/3 1/3 

1/3 1/3 1/3 




8 S 








S a S 7 S 8 











1/3 1/3 1/3 

1/3 1/3 1/3 

1/3 1/3 1/3 

1/3 1/3 1/3 . 





Line and 












Source 


Average Queue Sizes 








No. 


Utilization 


Si Si Sa S4 S-, Si 


Si Sa 


Ave. 


Max. 


1 


X = 0.239, 
P = 0.956 


9.7 10.1 13.8 6.7 15.6 14.6 


11.2 6.8 


11.1 


33-58 


2 


X = 0.229, 
P = 0.916 


5.6 5.6 6.1 3.0 5.0 6.8 


5.2 3.2 


5.1 


26-45 


3 


X = 0.213, 
p = 0.852 


2.4 2.6 1.9 1.1 2.4 2.3 


2.1 1.2 


2.0 


16-22 



The alternating priorities effect, due to domination, is demonstrated 
in Fig. 3. The total number of packets at the four queues of d (dotted 
line) and C 2 (solid line) were plotted against time. For the case p = 



2546 



THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1971 



^^^^^^^}Se^ ^wEl^ ^ 



X = 0.213 
p = 0.852 



.-7 100 



80 



60 



20 



SUM OF QUEUE SIZES AT C, = (S 1f S 2 , S 3 , S 4 ) 

SUM OF QUEUE SIZES AT C 2 = (S 5 ,S e , S 7 , S B ) 




50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 

TIME 



Fig. 3 — Simulated queue sizes. 



0.956 (high utilization) one can clearly see that only one dominating 
class is served at one time. When stations of C 2 take over control 
of the system the queues at stations of C i build up while at C 2 they 
are being depleted until the queue at S a reaches zero and C x can take 
over. f The average cycle time (time elapsing between two consecutive 
peaks of the dotted or solid lines) for this case was 110 time units. 
As the load on the system is decreased, the alternating priorities effect 
becomes less and less distinctive. For p = 0.916 (not graphed) the 
average cycle time reduces to 23 time units, and for p = 0.852 (see 
Fig. 3) alternations are very frequent and cycles are practically un- 
noticeable. 

Notwithstanding the complete symmetry within classes the average 
queue sizes at £ 4 and S 8 are consistently smaller than the average 
queue sizes at the other stations. The explanation of this phenomenon 
is as follows: At a moment when C, loses control to C 2 the queue size at 
St is zero while at $i , S 2 , and *S 3 it is greater than or equal to zero. 
From that moment on the queues of d build up at equal average rates 
until the moment control returns to C x (peak of the dotted line in Fig. 3). 
At that moment the queues start being depleted at equal average 

t Note that the queue size at Sb being zero is a necessary but not sufficient condition 
for losing control. 
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rates. Since the expected queue size at S t is smallest, it has a higher 
probability of being completely depleted first. This phenomenon tends 
to shorten the alternations cycle. 

The alternations cycle will be shorter when the number of rotating 
arms (loop time) is increased. Therefore large loops (local or regional 
loops) will be relatively more stable when high utilization occurs. 

The nature of the stochastic process used for generating packets 
at the stations of the loop is described in Ref. 5. 

An important aspect, not analyzed in this study, is the question 
of the order of stations in the loop. Clearly, the amount of traffic the 
loop can carry and the resulting congestion are strongly dependent 
on the specific order of the stations in the loop. In Example 1 we have 
assumed counterclockwise traffic direction. If we reverse the direction 
of traffic on the main line we shall get a different solution for the flows. 
The two solutions are compared in the following: 



Counterclockwise 


Clockwise Flow 


Flow Direction 


Direction 


x x = 0.4646 


x, = 0.5000 = X, 


x 2 = 0.4000 = X 2 


X a = 0.4000 = X 2 


x 3 = 0.5000 = X 3 


x a = 0.5000 = X a 


x 4 = 0.2333 


x t = 0.0917 


x., = 0.0508 


x 5 = 0.3667 



Total: 1.6487 1.8584 

Reversing the flow direction results in an increase of 13 percent in 
the amount of satisfied demand. 

In a practical situation not all orders are feasible. Still the number 
of feasible orders may be overwhelmingly large and an appropriate 
algorithm for determining best order is called for. An interesting pos- 
sibility is a double loop system where each station is connected to two 
loops with opposite traffic directions. This may increase reliability 
and enable better utilization by allocating traffic in an efficient manner. 
One possible allocation rule is shortest distance allocation where the 
loop to be used is the one with the shortest travel distance for each 
particular message (this is an example of a possible rule and is not 
proposed as an optimal rule) . 

The bounded linear complementarity problem presented by equa- 
tions (7) is of somewhat more general interest, bearing little relation 
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to the Pierce loop. In matrix form we have 

X + U - X, 

AX + Z m 1, 

C/ r Z = 0, U £ 0, Z ^ 0, Z ^ 0. 

Substituting Z = X — C/ we obtain the same set of relations in a slightly 
different form. 

AU-Z=A\-l = q 

U T Z = 0, £/" ^ 0, Z ^ 0, U £ X. (28) 

It is now possible to compare our problem to the Fundamental Problem 1 
treated by Lemke, 6 and Cottle and Dantzig. 2 

AU - Z = q 

U T Z = 0, U ^ 0, Z ^ 0. (29) 

The only basic difference between the two problems is that in our 
problem U is bounded from above while in Lemke's problem U is 
unconstrained. In this respect our problem is more general. The shift- 
and-load procedure is, however, fundamentally based on the specific 
structure of A and q and may not prove useful for a wider class of 
parameters. 
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APPENDIX 

In this appendix we outline a procedure for determining all dominat- 
ing submatrices of a given stochastic n X n matrix P. 

Step 1: Construct a matrix n = {t,,} where 

"l if p it > 0, 



*a = 



.0 otherwise. 



+ The term "Fundamental Problem" was coined by Cottle and Dantzig. 2 
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Step 2: Q = \qu\: = n. 

Do for j — 1, 2, • • • , n. 

Repeat until q ki , k = 1, 2, • • • , n, remain unchanged. 

ffw : = fo + S g« , * - 1,2, ••• ,n. 

N :a.-,-H 

(Note that all additions are Boolean.) 

If Si — > £,• then in the resulting matrix Q the element g, y = 1. 

Otherwise </,, = 0. 

Step 3: Construct a matrix Q <0> = {g-,} such that 

C = g!? J = m« ■ 

If S, does not belong to a communicating class then qlf = 0, 
j = I, 2, ■ ■ ■ , n. Otherwise gj? = 1 and *S, c C(i') if and only 
if g,-"' = 1. C(i) is closed (dominating) if and only if the z'th 
row of Q l0) is identical to the ith row of Q. 

REFERENCES 

1. Pierce, J. R., "Network for Block Switching of Data," unpublished work. 

2. Cottle, It. W., and Dantzig, G. B., "Complementary Pivot Theory of Mathemat- 

ical Programming," Mathematics of the Decision Sciences, Part 1, American 
Math. Roc, 196S, pp. 115-136. 

3. Kemeny, G. J., and Snell, J. L., Finite Markov Chains, New York: I). Van 

Nostrand Co., 1965. 

4. Avi-Itzhak. B., Maxwell, W. L., and Miller, L. W., "Queuing with Alternating 

Priorities," Opcr. Res., 13, No. 2, 1965, pp. 306-318. 

5. Hayes, J. F., and Sherman, I). N., "Traffic Considerations in a Pierce Data Loop, 

unpublished work. 

6. Lemke, C. E., "Bimatrix Equilibrium Points and Mathematical Programming, 

Mgt. Sci., 11, No. 7, 1965, pp. 681-689. 



